Atom loss from Bose-Einstein condensates due to Feshbach resonance 
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In recent experiments on Na Bose-Einstein condensates 
[S. Inouye et al, Nature 392, 151 (1998); J. Stenger et al, 
Phys. Rev. Lett. 82, 2422 (1999)], large loss rates were 
observed when a time- varying magnetic field was used to tune 
a molecular Feshbach resonance state near the state of pairs of 
atoms belonging to the condensate many-body wavefunction. 
A mechanism is offered here to account for the observed losses, 
based on the deactivation of the resonant molecular state by 
interaction with a third condensate atom, with a deactivation 
rate coefficient of magnitude ~ 10~ 10 cm 3 /s. 
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Experiments have been carried out recently Jl|,|| (see 
also the review S) in order to control the interatomic 
interaction underlying the properties of a Bose-Einstein 
condensate (BEC). One way to achieve this JTJJ^] is by 
tuning a magnetic field B to modify the two-atom scat- 
tering length as predicted for a B-dependent Feshbach 
resonance . The experiments carried out on a Na op- 
tical trap measured two distinct features: (a) the change 
in scattering length and large collisional atom losses with 
a slow sweep of B that stopped short of the resonant field 
Bq, and (b) a near-catastrophic loss of atom density with 
a fast sweep of B through the B region. Two groups of 
investigators |6]j7j have recently proposed a unimolecular 
mechanism to explain the latter feature as due to a fast 
sweep-induced transfer of population from the ground 
trap state of an atom pair to excited states. 

We present here another possible mechanism of atom 
loss, which is important when the field is changed slowly, 
and may play a role in the fast-sweep case as well. The 
temporary occupation of the vibrationally-excited molec- 
ular state m, coupled by the resonance to atom pairs in 
the condensate state a, is converted to a stable molecular 
dimer state d by a deactivating inelastic collision with a 
third atom: 

Na + Na^ Na 2 (m) (1) 
Na 2 (m) + Na -> Na 2 (d) + Na + AE (2) 

where AE is the excess kinetic energy released in the 
inelastic process. The deactivation states can be lower- 
lying rovibrational levels in the same spin state as the 
resonant level, or levels belonging to another spin state. 
While step (1), standing alone, is completely reversible, 
step (2) is irreversible, owing to the rapid ejection of 



the products from the condensate. The latter step usu- 
ally provides sufficient kinetic energy to freely escape the 
trap, and is characterized by a deactivation rate coeffi- 
cient 2j (in units of cm 3 /s). The two steps described 
above may be followed by a third one — secondary col- 
lisions of the products with condensate atoms. 

Except for 7, the magnitude of which can be extracted 
from fitting the experimental data, all expressions below 
can be derived analytically, using previously-determined 
parameters. A variational method was used to derive a 
set of three coupled Gross-Pitaevskii equations for the 
atomic condensate state tp a (r, t), the resonant molecular 
state ip m (r, t), and a representative deactivation state d. 
The latter can be decoupled out, leaving terms dependent 
on 7 in the remaining two equations (similar to those 
recently used by Timmermans et al. (3|^]): 

itUf>a= (j^P 2 + V * W + ^ B (*) + U a \<Pa?\ fa 

~^~U arn \ip m \ 2 Va+2g*(p* a (f m -ih^llfim^ifia (3) 
i%ip m = ^^P 2 + Vm (r) + /Ltm-B (t) + U m \<p m \ 2 ^J tp m 

|<£>a|Vm + g<p 2 a -ihj\tp a \ 2 <Pm (4) 

provided the deactivation state d never accrues a sig- 
nificant population. Here m is the mass of the Na atom; 
V a (r) and V m (r) are, respectively, the atomic and molec- 
ular optical trapping potentials (V m including the reso- 
nance detuning for zero magnetic field); r is the posi- 
tion in the trap; B (t) is the applied homogeneous mag- 
netic field; [i a and fj, m are, respectively, the atomic and 
molecular magnetic moments; and U a , U m , and U am 
are, respectively, zero-momentum atom-atom, molecule- 
molecule and atom-molecule interactions, proportional to 
the corresponding elastic scattering lengths. The cou- 
pling constant g responsible for the atom-Feshbach cou- 
pling is closely related to the parameter A that char- 
acterizes the strength of the resonance shape as a func- 
tion of the field B: 

\g\ 2 = 2Trh 2 \a a \[iA/m (5) 

where a a is the off-resonance triplet scattering length. 
Calculations |(^|| give a a — 3.4nm, /i = /i m — 2/i a = 
3.3^ B (where fi B = 9.27 x KT 24 J/T is the Bohr mag- 
neton), and A = O.OOlmT and O.lmT, respectively, for 
the two resonances observed at 85.3mT (853G) and 
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90.7mT (907G). These A values agree with the value 
measured for the 90.7mT resonance, and with the in- 
directly inferred order of magnitude for the 85.3mT res- 
onance. 

The analysis below neglects the kinetic energy terms 
in Eqs. (||) and (f|), in accord with the Thomas- Fermi 
approximation. This reduces (Q) and (Q) to a system of 
ordinary differential equations in tp a (r, t) and ip m (r, i) 
that depend parametrically on r. The resulting set of four 
real equations can then be solved numerically, for a given 
value of 7, starting with a Thomas- Fermi distribution or 
with a homogeneous (r-independent) distribution equal 
to the mean trap density. 

We have carried out such numerical solutions for both 
types of MIT experiments , using the homogeneous 
initial conditions. In the slow-sweep experiment the time- 
dependent magnetic field was linearly changed from an 
initial value of B to a final value closer to resonance, and 
then the density was measured by shutting off the trap 
field and letting the condensate expand. Repeating the 
experiment at various final values of the magnetic field, 
the loss was plotted as a function of this final value. It 
was noted that the results fit to a 3-body rate equation 
for the atomic density n(r,t) — \ip a (r,t) | 2 , of the type 
n w — K3T1 3 . 

Analysis of Eqs. (0, |]) shows that, out of the four real 
coupled equations for (p a and tp m , a single rate equation 
for the atomic density n(r,t) = \<p a (r,i) 2 can be ex- 
tracted, if the following condition (called here the "fast 
decay" approximation) 



\g\ 2 « h 2 j 2 n(r,t) 



(6) 



(i.e., n m (r,t) = \<p m (r, t)\ 2 <C n(r,t)) is observed. 
The resulting rate equation attains the nonlinear form 
n (r, t) = —F (n (r, t) , r, i), where 



F(n,r,t) 



6|<?|V 3 



[V{r)-^B{t)Y + {h in y- 



(7) 



and V (r) = 2V a (r) — V m (r). Eq. (0) can be expressed 
in terms of a collision cross section of the Breit-Wigner 
shape (dependent on the position r in the trap) for zero- 
momentum collisions, in which 2jn is the width of the 
output channel, and the width of the input channel is 
proportional to \g\ 2 . This observation establishes a link 
between the macroscopic approach used here and micro- 
scopic approaches, treating the loss rate as a collision 
process. 

It follows from Eq. (Q), that very close to resonance 
(where B (t) is within lfiT of a resonance), the behav- 
ior of Eq. (ffl) changes its form gradually from that of a 
3-body to that of an effective 1-body rate equation. How- 
ever, all values of B used for the slow-sweep experiments 
of the MIT group lie outside this narrow region, where 
the "off-resonance" condition 



holds. Given this condition, one can write Eq. (]?]) 
(to a very good approximation) in the 3-body form 
F (n, r, t) = K 3 (r, t) n 3 with 



K 3 (r,t) 



I2irh 2 \a a \-/A 
mil [B (t) - Bo (r) 



(9) 



where Bq(t) is the local resonance value of B. We shall 
refer to the combination of Eqs. @ and (||), leading to 
Eq. ([)]) , as the "3-body" approximation. This expression 
is similar to the one derived earlier |||| up to a factor 
3/2, which reflects the fact that 3 atoms are lost per one 
deactivating collision. 

For K 3 given by Eq. (^|), the rate equation may be 
solved analytically in the form 



n(r,t) = n(r,t ) 



1 + 24irh 2 \a a \Ajn 2 (r,t ) 



(t - t ) I (mfiB 2 tt 



1-1/2 



(10) 



h-yn(r,t) < \V (r) - 



(8) 



where B is the magnetic-field ramp speed, and the start- 
ing time to is measured from the extrapolated resonance- 
crossing time as t = 0. The graphs shown in Figs. 1 and 2 
pertain to the slow-sweep experiment conducted at MIT, 
and more particularly to the strong 90.7 mT resonance. 
This resonance has been approached from below with two 
ramp speeds, and from above — with one. The difference 
between Eqs. (||,Eo|) and the results of a direct numerical 
solution of Eqs. (|3|) and (Q) for all ramp speeds in the K3 
plots (Fig. 1) is so small, that the corresponding plots are 
indistinguishable. The remaining atomic density, at the 
moment t of stopping the ramp (see Fig. 2) , was also cal- 
culated, using the homogeneous-density initial condition, 
starting from a B value of 89.4mT on approach from be- 
low and 91.6mT from above, for which the initial mean 
densities were extracted from the experimental data [p|. 
The graphs clearly show a best fit with 7 of the order of 
10~ 10 cm 3 /s, which (given a density of about 10 15 cm~ 3 ) 
implies a deactivation time of ~ 10~ 5 s. 

The magnitude of the inelastic rate coefficient 27 
needed to explain the data is quite plausible. It is two or- 
ders of magnitude smaller than the upper limit allowed 
by the unitarity constraint on the S'-matrix. Accord- 
ing to this condition (hJ, in the limit of small momen- 
tum, 27 < hX/m, where A (the de Broglie wavelength) 
is in turn limited by the trap dimensions. This require- 
ment puts an upper bound of 2.5 x 10 _8 cm 3 /s to 27 in 
the case considered here. Furthermore, our estimate of 
10~ 10 cm 3 /s for 7 is consistent with calculations carried 
out recently |llj on the rate of vibrational deactivation 
in H2+He ultracold collisions, producing a comparable 
result in spite of the large energy gaps in H2 . 

It is interesting to compare the outcome of the deacti- 
vation mechanism described here with that of other loss 
mechanisms based on a unimolecular process (not involv- 
ing a third atom) [BUT]]; proposed for the second type of 
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MIT experiment with the fast sweep of B through reso- 
nance. To see this one can remove the imaginary term in 
Eq. (H|) , and replace in Eq. (g) with an n-independent 
r. In the fast-decay limit, in which n m <C n, both mech- 
anisms would yield similar results for K3 and n of the 
slow-sweep experiments, if similar decay times (~ 10 _5 s) 
were to be reckoned in both processes. 

The 3-body approximation does not hold very close to 
resonance, and is therefore inapplicable to a description 
of the fast-sweep experiment, in which the Zeeman shift 
was swept rapidly through the resonance, causing the 
most dramatic loss rates (see Rcfs. Nevertheless, 
maintaining only the fast decay approximation, Eq. (|^), 
a simple analytical expression can be obtained also for 
this kind of measurement. Assuming the magnetic field 
variation is long enough to allow taking the asymptotic 
time limits, reached under the condition /j,5B ^> hjn 
where SB is the range of variation of B on either side of 
the resonance, the ramp starting and stopping times can 
be extended to ±00. One then obtains at all positions r 



n (r, 00) 



n (r, 



1 + sn (r, —00) ' 



127T 2 % a | A 

m \B\ 



(11) 



This asymptotic result (as opposed to the one mentioned 
earlier with regard to the slow-sweep MIT experiment) 
is independent of 7. 

Assuming, as before, a homogeneous initial density 
within the trap, Eq. ([ll]) also describes the loss of the 
total population N (t) = J n (r, t) d 3 r. An analytical 
asymptotic expression for the total population can also 
be found when the homogeneous distribution is replaced 
by the Thomas- Fermi one (given in Ref. JT2I ) , in which 
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(12) 



where hq is the maximum initial density in the center of 
the trap. 

The analytical results of Eq. ( |ll|) , together with the 
direct numerical solutions of Eqs. (|3j) and (|4|), are com- 
pared in Fig. 4 with the results of the fast-sweep ex- 
periment pL^j- The numerical results show that if the 
deactivation time is not small enough compared to the 
sweep time, and the magnetic field variation is not large 
enough, the assumptions underlying Eqs. (^) and JT2| ) 
do not hold, and thus the condensate loss will depend 
on 7. These calculations show a rough agreement with 
the measured loss rates for values of 7 of the orders of 
10 _9 -10 _11 cm 3 /s, though these results are not as sensi- 
tive to 7 as were the fits to the slow-sweep experiments 
discussed above. 

The deactivation reaction Eq. (|2|) produces 4qnn m 
"hot" particles (atoms and molecules) in a unit volume 
per unit time. Traversing a distance & at a speed v, these 



particles create new hot particles in a cascading process. 
The density of hot particles may then be estimated as 
rih « 4jnn m b/v exp (ba (n + n m )), where a is the elas- 
tic collision cross-section, considered here identical for 
atom-atom and atom-molecule collisions. The rate of 
additional loss due to these secondary collisions is esti- 
mated as vanrih for the atomic condensate and v<jn m nh 
for the molecular condensate. This loss rate depends es- 
sentially on the product ab. A conservative estimate of 
this product is obtained by taking a = 0.5 x 10~ 12 cm 2 
(the elastic Na-Na cross section at 0.1K jig ], a typical 
vibrational deactivation energy), and b = 10~ 3 cm (the 
condensate-cloud radius). A more realistic value of this 
product may be higher on account of the larger axial 
condensate length, the rise of the elastic cross section at 
energies below ~ lOmK |b|, and a possible contribution 
of inelastic collisions. Two values of ab are therefore used 
in Fig. 4. This figure shows that the effect of secondary 
collisions, calculated for a homogeneous density, may ap- 
preciably increase the condensate loss, while the use of a 
Thomas- Fermi density [see Eq. (|l2|)] should decrease it. 

The authors are most grateful to Wolfgang Ketterle 
and Joern Stenger for helpful information regarding the 
MIT experiments. 
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FIG. 1. The 3-body rate coefficient (K3) vs. the stopping 
value of the magnetic field, calculated with 3 different values 
of the deactivation rate 7 (in units of cm 3 /s), on approaching 
the resonance from below or above. These are compared with 
the experimental results (Ref. Ramp speeds shown are 

in units of mT/s. 
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FIG. 2. The surviving mean density vs. the stopping value 
of the magnetic field, calculated with 3 different values of the 
ramp speed B (in mT/s). Other notations as in Fig. 1. 
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FIG. 3. Ratio of surviving trap population N to the ini- 
tial one iV for the 85.3mT (853G) and 90.7mT (907G) res- 
onances in the homogeneous-density approximation vs. sno 
(where the parameter s is defined by Eq. ( |l"l| ) and no is the 
initial density). The curves show calculations carried out for 
different magnitudes of the coefficient 7 (in units of cm 3 /s), 
and with the analytical result of the fast-decay approxima- 
tion Eq. ([H]) (dots). The plots for the 90.7mT resonance 
with 7 = 10 _9 cm 3 /s and for the 85.3mT resonance with 
7 = 10~ 10 cm 3 /s are practically indistinguishable. The calcu- 
lations are compared with the results of the MIT fast-sweep 
experiment (Ref. 0). 
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FIG. 4. 



Same as Fig. 
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3, but taking account of secondary 
collisions for 7 = 10~ a cm 3 /s and various values of the pa- 
rameter ab (in units of cm 3 ). The analytical results (without 
taking account of secondary collisions) for the Thomas-Fermi 
distribution Eq. jl^ ) (solid line) and the homogeneous distri- 
bution Eq. (hi]) (dots) are plotted as well. 
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